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^ : ABSTRACT 

It is shown that including form factors for the quark-pseudoscalar meson couphngs 
of the Georgi-Manohar chiral quark model allows one to obtain the leading off-shell 
dependence of vr" — 77 mixing (as predicted by chiral perturbation theory) from the 
effect of quark loops on the meson propagators. Implications for p° — uo mixing and 
for the effects on meson mixing contributions to few body charge symmetry breaking 
observables are also discussed. 



INTRODUCTION 

It is commonly held that the bulk of non-electromagnetic charge symmetry break- 
ing (CSB) in few-nucleon systems is due to isoscalar-isovector mixing in the interme- 
diate meson propagators of one-boson-exchange graphs. The dominant contributions 
of this type are associated with p° — uo and, to a lesser extent, n" — rj mixing [1-4]. Re- 
cently, however, it has been pointed out that some problems exist with the standard 
evaluations of these contributions. 

First, recall that the p° — uo mixing matrix element is obtained from experimental 
data on e^e^ -^ tt+tt^ in the p — uj interference region and therefore corresponds to 
q^ ~ m^. The standard analysis assumes this matrix element is g^-independent and, 
therefore, uses the experimental value unchanged in the NN scattering region, for 
which g^ < 0. Not only has the validity of this ansatz never been tested, but, in Ref. 
5, using a model in which p° — uo mixing is generated by an intermediate quark loop as 
a consequence of the inequality of up and down quark masses, Goldman, Henderson, 
and Thomas (GHT) raise the possibility of significant g^-dependence of the p° — oo 
mixing matrix element. 

Second, a recent evaluation [6] oi -k" — t] mixing to one loop in chiral perturbation 
theory (ChPT), shows that 1) the mixing matrix element is, indeed, g^-dependent 
(varying by ~ 20% over a range from g^ ~ rni to g^ ~ ~f^% comparable to that 
involved in the extrapolation of the p° — uo matrix element from g^ ~ m^ to the scat- 
tering region) and 2) the magnitude of the mixing, even without this g^-dependence, 
is less than obtained by the model analysis of Ref. 3c. The results of Ref. 6, being 
at one loop in ChPT, provide only the linear-in-g^ dependence oi ti° — r] mixing; the 
higher order dependence is unconstrained but will be small in the region of validity 
of the chiral expansion (|g^| ~ fni)- Since the convergence of the chiral expansion to 
one loop in this region has been extensively tested [7,8], and the framework of ChPT 
is a rigorous consequence of QCD, the results of Ref. 6 should be considered very 
reliable (up to an overall ~20% scale uncertainty associated with the treatment of 
electromagnetic corrections of pseudoscalar masses only to leading order in the chiral 



expansion [6,9]). 

The aim of the present paper is to use the constraints of Ref. 6 to test the type 
of model building which underlies the GHT approach. As we will see, the modelling 
does, indeed, succeed in reproducing the correct g^-dependence of 7r° — 77 mixing, 
which considerably strengthens the case of GHT that, as a consequence of the q^- 
dependence of the mixing, the dominant p° — uj mixing contribution to few-body CSB 
may be significantly different than previously thought. 

THE Q2.DEPENDENCE ol -n" - r] MIXING AND THE QUARK LOOP MODEL 

Let us begin with a few general remarks. These are necessitated by the, apparently 
widely held, view that mixing matrix elements of the sort we are discussing are 
expected to be "naturally" g^-independent. This is simply not the case. 

In fact the whole assumption underlying the meson exchange framework is that, 
at low energies, QCD reduces to an effective low energy theory involving only com- 
posite hadronic fields (meson, nucleons, deltas...). Let us assume, for the sake of the 
argument in this paper, that this assumption is essentially correct. Then, whatever 
the effective theory governing the interactions of these hadrons is, we know that, as an 
effective low energy theory, it will be described by a non-renormalizable Lagrangian 
in which all terms not explicitly forbidden by the symmetries of the underlying theory 
(QCD) occur. In particular, there will be terms in the meson sector of the theory 
which involve the quark mass matrix and higher powers of derivatives, which will 
naturally lead to g^-dependent meson mixing. The pseudoscalar sector, where we 
actually know something about the effective theory beyond leading order in the mo- 
mentum expansion, is one explicit example of this general principle. The standard 
ansatz, of taking, eg., the p° — uo matrix element to be independent of g^, is thus 
incompatible with the assumptions underlying the framework in which such mixing 
is to be used to generate CSB in few-body systems. Of course, without knowing the 
terms of effective Lagrangian involving the vector mesons beyond leading order, one 
does not know what the magnitude of the g^-dependence is, and hence whether the 



effect of the g^-running is important or not. Given the pseudoscalar result and the 
natural (QCD) scales involved, however, it is unlikely that the effect will be negligible. 
The fact that the vector mesons are much less point-like than the pseudscalar mesons 
further supports this contention. Q 

Let us now turn to the question at hand, namely, whether or not a quark-loop 
model of the GHT type is capable of reproducing the behavior of tt" — r^ mixing known 
from ChPT. To one-loop in ChPT, Ref. 6 shows that one obtains a g^-dependent 
vTs — TTg mixing angle given by 
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where tts, vrg are the unmixed states to which the physical 7r°, r] states reduce in the 
limit of isospin symmetry, Aqmo is the Gell-Mann-Okubo discrepancy 



A 



GMO 



{Am 



m,„ 



'iml)l{ml-ml), 



(2) 



^Note also that the non-renormahzable structure of the effective low energy Lagrangian will 
lead to a g^-dependence of the "renormalized" (in the sense of the quark mass/momentum 
expansion) self-energies for the composite fields. This means that the meson propagators in the 
effective theory may differ significantly from the form 

A(9) = z/(g2 - m^) 

away from the pole (i. e. especially in the spacelike region relevant to NN scattering in meson 
exchange models). The pseudoscalar mesons are special in this regard since, as one can readily 
see from the structure the effective Lagrangian [see e. g. Ref. 7], the g^-dependence of the 
renormalized self-energy enters first at two-loop order and hence will be very small, by the 
usual power counting arguments, in the region of validity of the chiral expansion. This is a 
special property associated with the chiral constraints of the pseudoscalar sector and will not be 
a property, e. g. of vector meson fields. 



/ is a parameter of the chiral expansion (equal to /^ in leading order) and 
{rnj^o — Tn\+)QCD is the contribution to the kaon mass splitting due to m^ ^ m^, 
obtained by correcting the observed splitting for electromagnetic effects. This correc- 
tion is usually performed using only leading order results in the electromagnetic chiral 
expansion (Dashen's theorem [10]). If this procedure were accurate, the corrections 
to (1), which enter only at 6th order in the low-energy expansion, would be small in 
the region of small g^, say |q'^| ~ rni. 

An old attempt by Socolow [11] to saturate the Cottingham formula for the kaons 
with K, K* states, however, indicates that there may be significant corrections to 
Dashen's theorem. This possibility is also discussed in Ref. 9, on the basis of the 
size of the relevant electromagnetic chiral logarithms. At present, both the kaon 
electromagnetic splitting obtained from Dashen's theorem and the larger Socolow 
value are compatible with other isospin breaking data [9]. The larger Socolow value 
would lead to a value of {jn\o — rn\-+)QCD 30% larger than that obtained using 
Dashen's theorem. There is thus an overall scale uncertainty in (1) (which uncertainty, 
however, can be reduced by improved experimental results, 
e. g. on isospin breaking in K^^). 

Our aim now is to see whether a quark-loop model of the GHT type succeeds 
in reproducing the constraints just discussed. Two assumptions underlie the GHT 
approach: first, that the behavior of the full low-energy meson theory can be obtained 
by integrating out quarks from a theory consisting of free mesons coupled to quarks 
and, second, that the effect of the quark loops can be adequately represented by 
keeping only those loops generated by the lowest order quark-meson couplings of the 
effective quark-meson theory using a) a monopole form factor at each quark-meson 
vertex, and b) free constituent quark propagators. The basic philosophy behind these 
assumptions is that the higher order terms in the effective meson Lagrangian arise 
from the non-pointlike, quark substructure of the mesons, and that incorporating 
this substructure in a way that (through the monopole form factor parameter. A) is 
capable of reflecting the actual meson size should allow one to reproduce the behavior 
associated with these terms. 



For the problem at hand we, therefore, start with a model which incorporates both 
quarks and pseudoscalar mesons. The model must, moreover, properly respect the 
approximate chiral symmetry with which the pseudoscalars are associated as pseudo- 
Goldstone bosons. A natural choice is the Georgi-Manohar chiral quark model [12,13]. 
The model describes an effective theory of massive constituent u, d, s quarks and the 
pseudoscalar octet, {tt"}, transforming non-linearly under chiral SU{3)l x SU{3)r 
according to 

e ^ e' = L^U+ = U^R+ (3a) 

q^q' = Uq (36) 

where ^ = exp{iX ■ 7r/2/) is the square root of the matrix S = exp{iX ■ vr//) , { A" } 
are the usual Gell-Mann matrices and / is as in Eqn. (1) above. S transforms linearly 
under SU{3)l x SU{3)r, 

S^S' = LSi?+. (4) 

In (3) and (4), U reduces to the usual flavor SU{3) matrix for an SU{3)y transforma- 
tion and is a non-linear function of L, R, {vr"} otherwise. In the chiral limit (in which 
the model is to be invariant under the transformations (3a), (3b)), if we consider only 
terms with zero or one derivative (the leading terms in the momentum expansion), 
the effective Lagrangian of the model is 

where the covariant derivative, D^, is given by 

D, = d,- tV, (6) 

and the vector and axial vector fields V^, A^ are defined by 
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In Eqn. (5), Qa is the constituent quark axial coupling, g^ — 0.75, and m is 
the hypothetical constituent quark mass, in the chiral limit. At this stage chiral 
symmetry is unbroken, the vr" are all massless, isospin is exact, and there is no tt^ — tts 
mixing. Chiral symmetry breaking is then incorporated using the standard techniques 
of effective chiral Lagrangians. The leading symmetry breaking term involving both 
quarks and pseudoscalars is 
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where /i is a mass scale related to the quark condensate, M is the current quark mass 
matrix, M = diag{m1,m'^,m1), A^ is a chiral symmetry breaking scale ~ IGeV and 
C is expected to be of order 1. C is in fact constrained by the observation that the 
terms zeroth order in vr" in Eqn. (8) produce a splitting of the constituent s mass 
from the constituent u, d masses 

(5mr = m^ - ^M T^ ~ 200Me\/ (9) 

X 

where the lowest order mass relations for the pseudoscalars have been used to set 
^iTLs — Tn\. For further details of the model, the reader is referred to Refs. 12,13. 

According to the GHT ansatz, what we are now supposed to do is 1) drop all 
additional terms in the chiral expansion for the quark, pseudoscalar, and quark- 
pseudoscalar sectors, and 2) hope to generate the effects of these terms by including 
a monopole form factor, F{k'^) = A"^ / {A"^ — k'^) , for each meson leg coupling to a quark 
line (with k the four-momentum flowing through the vertex on the quark line) and 
considering the effect of quark loops with free, constituent quark propagators. 

Such loops contribute to the pseudoscalar propagators in two ways. First, the 
axial vector coupling term in Eqn. (5) generates pseudovector quark-pseudoscalar 
meson couplings, which in turn generate two- vertex loops as in Fig. 1. (Because 
of the pseudovector nature of the couplings, these loops are proportional to q^ near 
g^ = and do not contribute to the pseudoscalar masses.) Second, the terms of Lb 
second order in the pseudoscalar fields produce tadpole diagrams, as in Fig. 2. These 
loop contributions are proportional to current quark masses and naturally lead to 



the usual relations between the pseudoscalar squared-masses and the current quark 
masses, 

m^+ = fiim"^ + ml) 
mj^o = n{ml + nQ 
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up to an overall factor which is A-dependent, as we will see in more detail below. 

Requiring that we reproduce the correct leading mass relations would, therefore, 
fix the value of A to be used in the GHT ansatz. Choosing A in this way then 
automatically also fixes the tts — vrg mixing angle to its correct, leading order value 
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The real test of the GHT ansatz is then whether or not it is able to reproduce the 
correct g^-dependence of 6*, as given in Eqn. (1). 

To evaluate ^(g^) in the model, we include the quark loops of Figs. 1,2. The tts — TTg 
inverse propagator then takes the form 
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In Eqn. (11), vtss , ttss , vrgg all include both tadpole and two-vertex-loop contribu- 
tions. The former are g^-independent, the latter g^-dependent (and beginning at 
0{q^) near g^ = 0). TTag, moreover, is proportional to (m^ — ml). For small g^ (the 
only values for which we have a constraint) we may write 
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where the ir^i are associated with the tadpoles and the tt^/ with the two- vertex loops. 
Then the angle which diagonalizes A^^, to 0{md — 'my) (the same order as the ChPT 
result in Eqn. (1)), is given by 
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where we have kept only terms up to 0{q^) in the expansion occuring on the second 
line of Eqn. (14). The appearance of the leading order expression as an overall factor 
in Eqn. (14) ensures that the scale of 6{q^) at g^ = is the leading order result, 
Eqn. (1), when A is chosen so as to give the leading order results for the pseudoscalar 
masses. There are, in fact, corrections to the leading order pseudoscalar mass expres- 
sions at one-loop in ChPT which would affect the choice of A, but these are not fixed 
numerically with great precision because of residual uncertainties in the higher order 
low-energy constants of Ref. 7. The resulting overall scale uncertainty in Eqn. (14) 
is not, however, a practical difficulty for three reasons. First, we know that the value 
should be close to that required to give the correct leading order mass formula and 
that this value of A gives the leading order result for 9, which will be close to the 
correct magnitude at g^ = even after one-loop corrections. Second, there is a range 
of uncertainty of 30% in the overall scale of the result Eqn. (1) as a consequence of 
the uncertainties in the electromagnetic corrections to the kaon splitting. And third, 
it turns out that the slope with q^ of the factor in brackets in Eqn. (14) is rather 
weakly dependent on A. We will, therefore, quote all results below with A fixed so 
as to give the leading order pseudoscalar mass relations. The potential error in the 
slope with q^ is ~ 10%. 

THE QUARK LOOP CONTRIBUTIONS 



To generate the loop contributions of Figs. 1,2 we require the relevant quark- 
pseudoscalar meson couplings. The pseudovector couplings required for Fig. 1 come 
from the third term in Eqn. (5) and are given by 
-9a /_ 7 A r,u 9a 



2/ 



(m7m75M - d-iiJilbd) d^T^s 



2VSf 



UJ^ll5U + d-ff,-f5d - 2s7^75s] S^tts 



a=3,8 \ J / 



(15) 



where the second line defines the couplings, gpy. The tadpoles of Fig. 2 are produced 
by terms which couple two pseudoscalars to a quark line. The contributions associ- 
ated with the vector current part of the second term in Eqn. (5) vanish by Lorentz 
invariance. The surviving tadpole contributions are then generated by those pieces 
of Lb (Eqn. (8)) second order in the {it"'}, the terms involving tts, vrg being 
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where the second line defines g^ad- 



From Eqn. (15) we find that the 2- vertex- loop contribution to the ab element of 
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the mass matrix, T^ioJpi'f)^ is given by 
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This integral is readily evaluated using standard Feynman parameter techniques. The 
resulting expression is somewhat lengthy and will not be quoted here since we require 
only the leading, O(g^), contribution. This may be obtained straightforwardly from 
the full expression, and is given by 
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Similarly, from Eqn. (16), we obtain the (g^-independent) tadpole contribution to 
the mass matrix, Tiah 
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(A2 - m2) (A2 - m2)2 
and, in the second line of Eqn. (19), we have retained only those terms linear in the 
current quark masses, m^. Eqn. (19), together with Eqn. (9), implies 
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which, as promised, are the correct leading order results providing A is chosen so that 

Q{A,m)=mlf/6mr- (22) 

RESULTS 

In order to employ the results of the last section in Eqn. (14), we require a value for 
the constituent quark mass, m. An extensive analysis of the meson [14] and baryon 
[15] sectors, which includes kinematic relativistic corrections, obtains m = 220MeV, 
5mf^ = 200MeV. As an alternate value for m, and in order to display the sensitivity 
to the parameters, we will also quote results for m = jItIt^o, which would put the 
threshold for the mixing matrix element in the correct location. Rewriting Eqn. (1) 
as 

QChPT^^2^ = ^(g2 = 0) [1 + cig2] (23) 

we have, from Eqn. (1), that Ci = 0.279^6^^^ 

In Table 1 we present the model results for ci. The errors quoted correspond to 
allowing A to vary over a range for which the coefficients of the leading order results 
in Eqn. (21) vary by ±20% about the value 1 (a typical variation associated with 
one-loop corrections in ChPT). As may be seen from the Table, this error amounts 
to ~10%. Ci also varies by ~5% between m = ^mj^o and m = 220MeV. A similar 
variation is produced by varying 5mf^ by ~10%. As we see from the Table, the 
agreement with the ChPT result is quite good. This gives us confidence that the GHT 
type of modelling may give reliable results for the off-shell dependence of mixing, at 
the very least to lowest non-trivial order in q^. 

It should be pointed out, in passing, that the GHT calculation only investigated 
the g^-dependence of the off-diagonal p — uo matrix element, TTp^;, in the inverse propa- 
gator, analogous to ttss of Eqn. (11). When one, however, extracts the p — a; element 
of the propagator itself and writes it in the form 
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(24) 



the g^-dependence of T^^J^^^{q^) arises not only from the proportionahty to vTpt^, but 
also from the g^-dependence of the p, u renormalized self-energies. The actual GHT 
predictions correspond to T^p^ and not vrJJ^^" and so should not be employed, in their 
present form, in few-body CSB calculations. 

The present calculation, however, suggests that the full version of the GHT calcu- 
lation for puo mixing should provide a reliable estimate of the off-shell behavior of the 
puo matrix element, at least in the region of small q^ . If one, further, wished to take 
the success of the calculation as evidence for the validity of the physical argument un- 
derlying the GHT ansatz, one might hope that the higher order g^-dependence would 
also be well-modelled. We have, of course, no way of demonstrating that this will be 
the case; however, from both the present results, those of ChPT, and from general 
principles, it is clear that a model which passes the test of satisfying at least one 
known constraint should be considered more reliable than one (the standard ansatz 
of a g^-independent matrix element) which fails it. This means that the few-body 
CSB contributions to such observables as the difference of n and p scattering lengths 
and the non-Coulombic A = 3 binding energy difference need re-evaluation. This can 
only be done reliably (i. e. in the regime where the modelling has been tested) to the 
extent that meson-mixing contributions are associated primarily with the region of 
g^ for which the leading g^-dependence of the mixing matrix element is dominant. 
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TABLE 1. The slope of 0{q'^) with g^ as a function of m, 5m'l 



m(MeV) 5m^°"(MeV) ci{GeV~'^) 

202.5 175 .223 ± .026 

200 .234 ± .026 

225 .244 ± .029 

220 175 .211 ± .024 

200 .222 ± .025 

225 .231 ± .026 
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Figure Captions 

1. Pseudovector-coupling- induced contributions to tlie meson propagator matrix. 

2. Tadpole contributions to the meson propagator matrix. 
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